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A Computation of Cyclic Cubic Units ' 

Harvey Colin and Saul Gorn 2 

The paper is a report on a tabulation of units performed on the EDVAC at the U. S. 
Army Proving Grounds in Aberdeen, Maryland. The algebraic number fields involved 
were 45 cyclic cubic fields of discriminant I 2 , where I is one of the primes of form 3ra+l 
from 7 through 499. The object of the search was the discovery of units through a specific 
method, an algorithm based on Minkowski's geometric number theory but particularly 
amenable to a stored-program computer. In the computation, combined use was made of 
integral arithmetic and decimal arithmetic but with careful error analysis. 

1. Introduction 

An essential part of the lure of algebraic number theory consists of tabulations of almost 
unpredictable irregularity. The advanced theory has somehow surpassed the scope of numeri- 
cal data without perceptibly improving the power to predict such data. As a result, for a 
small calculation an algebraic number theorist would not forego the personal satisfaction that 
results from the many ingenious, fortuitous, and deeply meaningful combinations of integers 
that he would encounter. Yet to really appreciate these vital irregularities he needs longer 
tabulations, with the inevitable requirements of uniformity and efficiency. Here the modern 
electronic digital computer can be used profitably. 

The discussion of the problem that follows will be primarily from the point of view of the 
computer program rather than from its theoretical development. Thus, for instance, the field 
will have to be regarded not as a set of elements satisfying certain axioms, etc., in the manner 
of Dedekind and Steinitz, but as an algebra with a definite basis and definite "structure con- 
stants" for multiplications. As a matter of historical fact the exhibiting of a basis was gen- 
erally far from trivial and lagged considerably behind existence theorems [1, 10} 3 

2. Description of Fields 

For any prime /( = 3ra+ 1) a cyclic cubic field will be defined. Our interest centers around 
algebraic integers of such a field, written as follows (owing to the presence of the so-called normal 
basis) : 

£=aa>i + &<o 2 +ca} 3 — (a,b,c). (1) 



Here a, b, c are rational integers (or coordinates), and co ly co 2 , co 3 (or the normal basis) are 
three roots of the defining (normal) equation [8]. 


the 


F(w) = u 3 + u 2 -ma-n = 0. 


(2) 


Here n is defined as follows: First, 4/ is decomposed (uniquely) into the combination 




4l=A 2 +27g 2 , 


(3) 


where the signs of A and g are specified by the conditions (in integral k) 




A=3k+1, #>0. 


(4a) 

The 
atory, 

round 
ed the 


1 The computing time was provided under contract DA-028-ORD-12332 with the United States Army Office of Ordnance Research 
first named author was then at Wayne University, Detroit, Mich., and the second named author was at the Ballistics Research Laboi 
Aberdeen, Md. The first named author is currently at Washington University, St. Louis, Mo. 

The preliminary report (reference [5]), was cleared for publication by the Ballistics Research Laboratory on April 26, 1956. The backg 
description, the illustrative examples, and the final tabulation were prepared at the National Bureau of Standards. This work has form 
basis for a further calculation currently on the SEAC. 

2 Now at Moore School of Electrical Engineering, Philadelphia, Pa. 

3 Italicized figures in brackets indicate the literature references at the end of this paper. 


4.">2248— 57 1 155 





Then, finally an integer n always exists, for which 

27n= (A+3)l-l. 



(4b) 



All algebraic integers under discussion are real. 

The most general basis of the field would be given by an integral unimodular combination 
ti= y ^a ij u j ij = 1,2,3, det|a^| = ±l. There is essentially no other basis [9]. 

The particular designation of roots is standardized as 



Thus we can distinguish among conjugates £ (given above), and £', £" given by 

£'=c<oi + aco 2 +6w 3 = (c,a,b) 
^ f = bo)i J rcoo 2 + ao) 3 = (b,c,a) 



(5) 

(6a) 
(6b) 



(formed by cyclic permutation of co x , « 2 , w 3 ), so that £'"= £, etc. As a matter of notation, the 
subscripts denote conjugates only in reference to the basis co*. Later on £ f shall denote (unre- 
lated) integers of the field. 

It was mentioned earlier that the discriminant of the field or the determinant 



(7) 



has the value I 2 . This fact (which can be verified directly from eq (1)) is not used directly 
in the machine computation but uniquely identifies the field in question [6], 

The multiplication of any two numbers in the field can be expressed in integral arithmetic; 
e. g., 

£i=(ai,&i,Ci), £2=(a 2 ,6 2 ,c 2 ), (8) 

then one can find the explicit representation of £j£ 2 in terms of the normal basis by means of 
the structure constants a u ft, y u a 2f ft, 7 2 , defined by 



0>i 


co 2 


0?3 


W 2 


C0 3 


COi 


0> 3 


Wi 


co 2 



and 



o>i = a 1 co 1 + ft co 2 +7ico 3 



co 2 w 3 = a: 2 Wi + ft w 2 + 72^3- 



The explicit formula is, then, e. g. 
where 



A = aiaia 2 + 716162 + ft^ 2 + ft (a!& 2 + a 2 60 + a 2 (&iC 2 + 6 2 c0 + 72 (ci^ + c 2 aO , 
5=fta 1 a 2 +a 1 6 1 6 2 +7 1 CiC 2 +7 2 (a 1 6 2 + a 2 6i) + ft(6iC 2 +6 2 c 1 ) + a 2 (cia 2 + 
(7= 7i^i«2 + ft 6i6 2 + odCiC 2 + a 2 (a x b 2 +0260 + 72 (61C2 + 6 2 c0 + ft {c x a 2 + c 2 a x ) . ^ 
This product can therefore be computed in the following tabular form: 



(9a) 



(9b) 



(1^0 


a x 


ft 


7i 


6162 


7i 


«i 


ft 


CiC 2 


ft 


7i 


«i 


&i6 2 +a 2 &i 


ft 


72 


«2 


6iC 2 +6 2 Ci 


QJ 2 


ft 


72 


Cia 2 -\-aiC 2 


72 


«2 


ft 




(4 


£ 
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These structure constants are given explicitly by the formulas: 

m + k+1 g+m—a 2 —g + m — a 2 

OL 2 = o "' ft = o > 72 — ' 



■1 



(9c) 



3 2 ll 2 

. g—m — a 2 —g—m—a 2 
a 1 =a 2 —m + lt ^i=- g > 7i= 9 J 

They are listed in table 2 (see p. 166). 

In particular, the norm of £ is defined for (1) by the integer 

^£)=££Y'=-^(a+&+c) 3 -^(^ (10a) 

and the trace is defined by 

*({)=$+{'+{"=- (a+6+c), (10b) 

so that the defining equation of £, £', or £" is 

with the further designation 

i(£) = t(tt') = -m(a+b+c) 2 +l2ab. (10c) 

Note that if a+6+c=0, then iV(£) is divisible by I. Thus, applying this information to 

{-{'=((a-c),(6-c),(c-a)), 

we find iV(5— £')=(£— £')(£' — £")(£" — £) ^ an integer divisible by I. 

3. Units 

Units are algebraic integers of norm equal to ±1. They are written as 

r)=U 1 Wi+U 2 a>2+U z a)z=(Ui,U 2 ,U3) ) (11) 

where, as before, u u u 2 , u z are integers. The most general unit has the form i^i^ 6 , where a 
and b are (positive, negative, or zero) integers, whereas 770 is a specially designated fundamental 
unit [6, p. 19]. The only fundamental units are ±rj^ 1 , ±(r}' ) ±l , ±(Vo') ±1 - In what follows, 
units will be normalized (by a change of sign) to have, conveniently, only norm +1. 

If we take eq (10a) and set iV(£) = l and (a,b,c) = (u u u 2 ,u 3 ), we obtain a cubic equation in 
these three latter unknowns, which we are in effect solving. Now such an equation (unlike, 
say, the Pell equation), has very little intrinsic interest. The interest in units is derived en- 
tirely from their role in unique factorization . 

To give an example, when £=19, one finds (as a byproduct of the main calculation) that 

NM=N(2a> l -a> 2 )=7. (12) 

This does not mean that 7 has two different factorizations, although seemingly 

7 = coiCo 2 a> 3 = (2 CO] — co 2 ) (2o> 2 — co 3 ) (2co 3 — coi). (12a) 

The reason, of course, is that the factors can pair off as associates, or 

(2co 2 — co 3 )/coi= (— o>i — 2co 2 ), (12b) 

a unit of the field. To see this, in (12a) calling one set of factors £, £j, £i' and the other 
£2, £2, £2, respectively, we must verify that one of the ratios £ 2 /£i> £2/^1, (a/fi ^ a unit. Thus 
for example, in accordance with the rules of section 2, we compute 
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and find that A*, B*, and C 1 * are now divisible by iV(£i). In practice, the whole computation 
would be performed modulo N(£ l )( = 7 here), using the structure constants of table 2. 

Thus the preliminary step in understanding factorization in fields becomes the recognition 
of units. Traditionally, one creates, by trial and error, combinations such as (12a), and one 
hunts for units, usually with better luck than one can completely explain [4]. 

The method used here will be very systematic. It is an algorithm for creating a surfeit 
of algebraic integers with small norms (including norm unity). In the process, a unit, indeed 
a fundamental unit, is inevitably produced. (Unfortunately, machine limitations forced the 
weakening of the method, as is explained in section 6.) 

4. The Sign-Discrimination Algorithm 

One operates with so-called reduced 3X4 matrices of algebraic integers [2, section 4]: 



r h 



<t>= 



b & ^ 



& & & & 

t" t" t" t" 

The columns represent conjugates of algebraic integers of sum zero 



(14a) 



(14b) 



and such that any three of the four & form a basis. The further, and most vital, condition is 
that the sign pattern be 



sgn fa 



■+ 



- + 
+ +J 



(14c) 



either as <j> now stands or under some rearrangement of columns. 

The algorithm is a method of generating such matrices in chains. A matrix fa is said to 
be a neighbor of <j> if it is formed by adding one column of <j> to another column of 4> and at the 
same time subtracting it from a third column of <£ (the three columns being different). Thus 
there are 4! possibilities to consider corresponding to the distinct triples (jij 2 jz), where 1 <j<4. 
For instance, the operation 



*[341fa 



(15a) 



can be understood to mean that the " third column of <j> is increased by the fourth column, which 
in turn diminishes the first," or, writing just the first rows, 



01= (£i— £4, £2, £3+ &, £4). 



(15b) 



(The relation (15a) can be read backward as an operation on fa.) Of the 4! possible neighbors, 
either 3, 4, 5, or 6 will be reduced, depending on inequalities satisfied by the elements [2, section 
10]. 

We next define the conjugate </>' of a matrix <f> to be the matrix formed by replacing the 
first row £ z - by its conjugate £* and forming the remaining rows from (■{ instead of £*. The 
conjugate <j>' will not be reduced as it stands, but a rearrangement of the columns (in this case 
the interchange of first and second) will render it reduced. Likewise, we define 4> n '. Note 
4/"=<£ again. Finally, two matrices fa and fa are said to be proportional to one another if, 
for </>* = </>!, <j)[, or <£'/, and numbers 77, ?/, 7?", 
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<t>2 = 



Srj 

7}' 

10 V'J 



0i* • (16) 



It is easily seen that 77, ?/, nf ] n are necessarily the conjugates of a unit, possibly 1, having 
norm +1. (In some cases <j> x is proportional to itself, leading to interesting possibilities related 
to ramified primes but beyond the scope of the present discussion [2, section 22].) In all cases 3 
proportional matrices will have the same set of norms (possibly under rearrangement). 

The algorithm operates as follows: The initial matrix is taken for convenience (writing 
only the first row) as 

0o=(«ij w 3, »2, 1), n>0, (17a) 



with norms respectively 



or- as 



with norms respectively 



(n,n, n, 1), (17b) 

0o=(<*>i + a>2, «8+«i, — coi, 1) n<0 (17c) 

(— n+m, —n+m, —n } 1). (17cl) 



Branching out from </> , one forms all reduced neighbors of matrices present, terminating a 
branch when the matrix is proportional to one already present. Indeed, the number of reduced 
matrices is finite to within proportionalities [2, section 8]. The factors of proportionality 
77 are units among which can always be found the fundamental unit [2, section 18]. 

5. Modified Algorithm 

The original algorithm, as just described, is a branching algorithm in which each </> may 
have more than one successor, and each c/> is compared with all those preceding. Because 
each </> contains 12 components and there could be more than 50 produced, the internal storage 
of the EDVAC would be taxed before the algorithm had gone very far. The branching 
algorithm was therefore modified to a form in which only one successor is chosen for each <j>, 
and the comparison (16) is always between <£ and the most recent </>. 

The modification consisted in trying only one-third of the 4! possible neighbors, namely, 

[123], [132], [423], [432], [214], [241], [314], [341] (18) 

characterized by the fact that the second and third indices are either 2,3 or 1,4 in some order. 

Either one or two neighbors of this type occur, but the EDVAC takes only the first that occurs 

in the program, thus we have a chain-type algorithm rather than a branching algorithm, with 

4 
the further property that Sl&l decreases at each step [2, section 17], whereas the corresponding 

1 

sums for the second and third row (conjugates) increase or remain the same. This assures 
us that when the inevitable proportionality occurs, the unit will not be 1 (but it need not be 
a fundamental unit either). 

As an illustration (see table 1) 4 the matrix sequence <I> 1 is reproduced for 1= 19. It is 

<M241] ^[132] 4> 2 [314] 4> 3 [314] <fc[314] <fe[423] 4> 6 [341] </> 7 [123] tf> 8 [123] </> 9 [123] tf> 10 [432] *n[214] </> 12 . 

The EDVAC discovered that </> is proportional to (some rearrangement of) </> )2 . The specific 
numerical values are of some interest. We therefore calculate 



i A full explanation of table 1 is given in section 7. 
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o)!= 2.507 018 643 
W2= — 1.221 876 162 

co 3 =-2.285 142 481, 

which (except for minor modifications in section 7) is the accuracy presented by the EDVAC. 
Thus 

<t>o= ((1,0,0) (0,1,0) (0,0,1) (-1,-1,-1)) 

or 

/ 2.507 018 643 -2.285 142 481 -1.221 876 162 1 \ 

<£<>=( -1.221 876 162 2.507 018 643 -2.285 142 481 1 j 

\-2.285 142 481 -1.221 876 162 2.507 018 643 1 / 

<t>u= ((7,5,5) (4-3,6) (2,6,-1) (-13-8-10)) 

or 

/ 0.014 037 286 -0.017 151 828 -0.032 077 205 0.035 191 747 \ 

</> 12 = { -7.443 752 324 17.010 034 653 -18.661 625 853 9.095 343 524 J 

\— 9.570 284 962 -23.992 882 825 11.693 703 058 21.869 464 729 / 

The EDVAC then "recognized" that the first column would serve as a proportionality factor, 
i. e., if we rearrange the columns of X2 to form 

*?2=((-13,-8,-10), (2,6,-1), (4,-3,6), (7,5,5)), 
then 



*?2 = 



V 
.0 i" 



0u 



where 77 = 0.014 037 286, V= — 7.443 752 34, and t?"=-9.570 284 962. In this case the 
decimal accuracy would be easily sufficient for "recognition," but integral arithmetic was used 
(see section 8). 

The unit produced, t]— (7,5,5), is not fundamental, but, early in the course of the compu- 
tation, for </>!, £ 2 happened to come out as rj = ( — 1,-1,0), a fundamental unit (see section 10). 
Actually, 

This can be easily verified from the multiplication scheme (9), e. g., (rjo) 2== (— 6,— 3, — 5), 
(r/o)- 1 ^ (2,2,1). 

6. Numerical Accuracy 

As the algorithm was first conceived and applied in experiment [2, section 11], decimal 
accuracy seemed sufficient because the symmetric functions could be correctly evaluated by 
the nearest integer when necessary. In a long computation, however, the algorithm can go 
astray in many ways: 

(a) The algorithm may incorrectly discriminate the reduced neighbors from the others. 
This is bound to happen because the subset of the reduced neighbors chosen by the modified 

4 

algorithm just mentioned has S |£*l approaching zero monotonically. The error in the quanti- 



ties £ t can soon be bigger than the quantities themselves, making a simple discrimination in 
sign untrustworthy. 

(b) It may make uncertain the decision as to whether <t> and <£* are proportional; for 
example, it may have to be decided that (using an obvious notation), £i/£i = ^2/&=kl&=^l^ly 
in which all the quantities involved are in error. 
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(c) It may of course make the approximations to the all-important resulting units untrust- 
worthy; for these are the ratios discussed in (b). 

A vital and desirable modification, to control the accuracy of the computations, consisted 
in carrying out all arithmetic operations exactly by using integral arithmetic; (explicit formulas 
for the integral arithmetic operations on two such integers are given in section 2). Furthermore, 
all quotients were removed from the algorithm (see section 8). 

It would also have been quite easy in principle to have handled the sign discriminations 
in integral arithmetic. For instance, to tell if two algebraic numbers £ x and £ 2 agree in sign, 
all we need have done is test to see if the three conjugates ^ 2} £i&, and £i£ 2 ' are each positive 
(total positiveness), which, by Descartes' law of signs, is a matter of seeing that the defining 
eq (10) for (£i£ 2 ) has (integral) coefficients with alternating sign. At the time, however, the 
EDVAC internal storage was not sufficient to handle the enormous integers, so that an auto- 
matic error analysis was used instead. This was the only departure from integral arithmetic. 

7. Input-Output and Error Controls 

With the exception of the following operations, the problem was run completely internally, 
with no human intervention: 

(a) For each /, the structure constants in table 2 were computed in advance by hand 
(merely to save space in the memory) through the a 2 column. (The entire table is given for 
convenience of hand checking.) 

(b) For each I, the zeros (and maximum errors) of the polynomial (2) were computed in 
advance by machine. The polynomial was scaled by dividing its zeros by m and carrying 
through the computation in fixed point. The root-finding method for the scaled polynomial, 

fM^-)'+^(-)'—(-)-% (19) 

is the bisection method (Horner's method in binary), for appropriate isolating intervals. The 
maximum error in the coefficients of the scaled polynomial, for EDVAC, is 2 -44 , whence the 
maximum error in evaluating the scaled polynomial is 7X2 -44 ; if, then, we use e=2~ 41 , the 
machine will provide an automatic error analysis in solving for co* by finding the zeros rf of 

/(co) + 2e, r z ~ of /(co) — 2e, and setting 

Wi rt+rl rt—n 



m 2 2 

Because the numerical values of the zeros co z are only used in the main algorithm to discriminate 
the signs of linear homogeneous expressions, aaj 1 + 6co 2 + cco 3 , the scaled values Wj/ra and their 
error bounds e t were used directly. 

(c) The input and output of the algorithm itself was by means of IBM cards, in decimal. 
(This required internal conversion on input and output because the EDVAC arithmetic is 
binary.) For each I two IBM cards were read in. One contained the integers /, m, g, a 2 , n. 
The other contained the integer I together with the (10-digit) decimals wi/m, (a 2 /m, w 3 /m, and 
the errors in «i, co 2 , o> 3 . 

The output cards were as follows: First, there were copies of the two input cards; 
second, for each </> u as soon as it was produced, four cards were punched with 7 = 1,2,3,4, 
respectively, indicating the column of <j> t (prior to possible rearrangement) being described. 
The algebraic integer ^ was described by the integers /, i, j, a h h h c h N(£), t(£), 7(f), where N, 
t, 7 refer to the defining eq (10) and 

Third, there would be a card containing the integers (see (11)) /, i,j( = 5), u u u 2 , u h iV(f ) ( = 1 ) , 
t(£), 7(f) if the algorithm runs to completion to produce a 4> t proportional to O with factor 7?. 
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Table 1. Input-output cards 
The material in this table is presented in the form in which it is printed out by EDVAC. 
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OUTPUT CARDS 
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21 
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Table 1. Input- output cards — Continued 



output cards — continued 
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Trie algorithm is then applied to the next case automatically. A complete set of input-output 
cards is listed in table 1 for Z=19. 

On the other hand, error controls could stop a case (or the current value of I) before it has 
run to completion in the following ways: 

(a) The first error control computes the numerical approximations to fr, £ 2 , £3, £4 from their 
components a h bj, Cj and values of «*. Thus the maximum error in ^=^a> I +6^co 2 +^w 3 is 
€ = |a i |€ i + \bj\e 2 + \cj\ez, so that the (scaled) condition €i/m<e would mark the sign discrimination 
on £, as untrustworthy. In this case an extra card was printed containing /, i, 5, £*, £*> £*, £% with 
£f being the numerical approximation. The computation proceeded to the next I automatically. 
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(b) The second error control checked each computed norm to stop the machine if any neg- 
ative norms appeared or if any </>, failed to have one of the eight neighbors (18). (None did.) 

(c) The third error control checked the imminence of overflow by testing to see if the 
largest term appearing in the norm calculation, n(a-\-b+c) 3 , becomes too large (>2 44 ) to 
fit in the registers. In case of this eventuality, the machine proceeded to the next case, 
printing no signal card. These controls required very little storage. 

8. Neighbor Formation and Comparison 

What remains is now to describe the portion of the program that is strictly internal, and 
which, by itself, would not come to the attention of the machine operator. The two main 
steps of the algorithm, are neighbor formation and comparison. 

The neighbor formation consists of the formation of </> z+ i from 4> t and, ultimately, the replace- 
ment of <\>i by <t>i+\. Here, essentially, each neighbor operation of type (18) was tried in succes- 
sion on <j>i, and the sign discriminations were made to test if the pattern (14c) were valid. 
Here the numerical estimates were kept in double precision, the integral and fractional parts 
being stored separately (each consisting of at most about 13 decimal digits). 

The comparison of <l> and </>*■ is more complex. First of all, the cases 1=7 and 13 were 
excluded by virtue of the fact that \n\ = l (so that + coi or — wi is a unit anyway). In the other 
cases (Z^19), the set of norms in wo (see eq (17b, d)) contains exactly one norm of value 1. 
Thus a <t>i is not conceivably proportional to fa unless some column in 4> t is a unit (say) rj. A 
further necessary condition is that all the norms in <t> t match those of <£ G - When this was satisfied 
the machine multiplied the top row of # by 77, V ', and 77", in turn (using (9a,b)) testing at each 
stage to see if <fo or a conjugate was obtained. (In this last step, the possible column rearrange- 
ments are fewer if there is only one norm in <£ equal to 1.) 

It is seen that the machine part of the algorithm, at no time involved division (even such 
minor calculations as division by 2 were performed by shifting the binary numbers). 

9. Supplementary Hand Computations 

The restriction that the whole process be contained in a memory of 1,024 cells proved to 
be very great. As a result, in the EDA r AC program provision was made only for double 
precision numerical estimates (the integral part and the fractional part being stored separately) 
and single precision integers. Thus only a few cases ran to completion, but only 18 cases 
failed to produce at least one unit in some <f> t (with or without giving rise to proportionality). 
In these cases a hand computation produced units from the algebraic numbers automatically 
computed. The general scheme for these hand computations was the following: Suppose the 
machine has produced two nonconjugate num.bers £ x , £ 2 with same norm, preferably a small 
prime p. If p is uniquely factorable in the field, then & must be associated with a conjugate 
of £ 2 , as we can test in the manner described earlier in section 3, completely in integral arith- 
metic. If the norm is not prime, more complicated combinations must be tried, but they are 
no different in principle. 

At least half the hand-computed units were obtained quickly in this fashion. The cases 
that did not yield to this method directly did so when more than two integers were used in 
proper combination. 

10. Fundamental Units 
In each case those units r\ were noted for which 

max t(ri),t(ri) (20) 

over all units produced in the algorithm was a minimum.. In accordance with arguments given 
by Hasse [6], this would be the fundamental unit if one were present in the list. 

By a remarkable coincidence, almost simultaneously with the completion of the runs 
(Aug. 1954) Peter Swinnerton-Dyer kindly communicated an independent table giving the 
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traces t, 1 of the fundamental units for 7<Z<300. A comparison of the tables showed perfect 
agreement except in the three cases, Z=373, 379, 463. In just these cases the (hand) computed 
units had positive traces for t and 7, which should have suggested that the units computed 
from the ED VAC output are squares. 

In the listing in table 2, the components and traces of both a unit and its reciprocal are 
shown in each case (here t{r})=t(l/rj)). In those cases where the units were computed by 
hand, (*) is shown. In general, regardless of how the units were produced or derived, the 
symbol "f" indicates the fundamental unit (and reciprocal) and "ff" is reserved for cases where 
both fundamental unit and reciprocal were produced by machine. In the three exceptional 
cases, the units, denoted by "s", are squares of fundamental units. 

In order, however, to give some kind of reasonable impression of the amount of labor, 
particularly in view of 13-digit coordinates (!), we should point out that the number of trials 
in testing fundamental units is easily 

0{log[max^),/,(^)]/log/} 2 . (21) 

To see this very quickly, note that if rj is not a fundamental unit, we can write (recalling 
vvv =1), 

v'=Wonvo) b vo b (vo) a - b y> 



--Vo b (vo) a - b [, 



where r? is some (as yet unknown) fundamental unit, and a and b are rational integers. By 
taking reciprocals, or conjugates if necessary, we can assure ourselves that a>b>(). Then 

ri b W) a =Vo M , (22a) 

where 

M=a 2 -ab + b 2 . (22b) 

Letting max (\rj\, |r?'|, |*?"|) = H, max (|iy |, |^|, |^J|) = Ho, we find 

H a+& >H M , (23a) 

or since (a+b)/M<2a/M<2a/(4a 2 /3)=Z/(2a), 

H <H (3/2a) . (23b) 

But, by the concluding remarks of section 2, 

1<W)— ^oll^o — Vo I ho — *7o|<8H 3 
and 

H 2 <„ 2 +(V) 2 +0?") 2 =*M 2 -2fa); 



thus 



-2 log H ~ log [//8] 



This verifies the order of magnitude (21). 5 

The problem of showing the units rj to be fundamental now resolves itself to checking on 
whether any of the units r) b (r)') a is a perfect M power, where a and b come from the finite set (24). 
But this check is somewhat laborious, even to program for the EDVAC, and owing to the 
convenience of Swinnerton-Dyer's table, it was never carried out. 



5 Much stronger inequalities are given by Hasse [6] using more specialized notation and information. 

165 



Table 2. Structure constants and units 






Prime 


Structure Constants 


I 


m 


k 


g 


n 


«2 


ft 


72 


ai 


ft 


Yi 


7 


2 





l 


1 


1 


1 





-2 


-1 


-2 


13 


4 


-2 


l 


-1 


1 


2 


1 


-4 


-2 


-3 


19 


6 


2 


l 


7 


3 


2 


1 


-4 


-4 


-5 


31 


10 


1 


2 


8 


4 


4 


2 


-7 


-6 


-8 


37 


12 


-4 


1 


-11 


3 


5 


4 


-10 


-7 


-8 


43 


14 


-3 


2 


-8 


4 


6 


4 


-11 


-8 


-10 


61 


20 





3 


9 


7 


8 


5 


-14 


-12 


-15 


67 


22 


-2 


3 


-5 


7 


9 


6 


-16 


-13 


-16 


73 


24 


2 


3 


27 


9 


9 


6 


-16 


-15 


-18 


79 


26 


-6 


1 


-41 


7 


10 


9 


-20 


-16 


-17 


97 


32 


6 


1 


79 


13 


10 


9 


-20 


-22 


-23 


103 


34 


4 


3 


61 


13 


12 


9 


-22 


-22 


-25 


109 


36 


-1 


4 


4 


12 


14 


10 


-25 


-22 


-26 


127 


42 


-7 


2 


-80 


12 


16 


14 


-31 


-26 


-28 


139 


46 


-8 


1 


-103 


13 


17 


16 


-34 


-29 


-30 


151 


50 


6 


3 


123 


19 


17 


14 


-32 


-33 


-36 


157 


52 


-5 


4 


-64 


16 


20 


16 


-37 


-32 


-36 


163 


54 


8 


1 


169 


21 


17 


16 


-34 


-37 


-38 


181 


60 


2 


5 


67 


21 


22 


17 


-40 


-38 


-43 


193 


64 


-8 


3 


-143 


19 


24 


21 


-46 


-40 


-43 



Result of algorithm 



No. of 

matrices 



special 
special 

12 comp 
30 comp 

13 comp 
32 over 
32 over 
19 comp 
34 over 
39 over 
50 over 
46 disc 
37 disc 
43 over 
36 over 

45 disc 
36 over 

46 disc 
22 over 
43 disc 



No. of 
units 



special 
special 

7ff 

3f 

3ff 

If 

2f 

If 

2f 

9f 
13f 

3f 

If 

If 

8f 

If 

If* 

8f 

If 

If 



Components of unit produced: 

71 ~ l — U\03i + W2CO2 + UzOiZ 



u x 



1 

1 

-1 

1 

-1 

2 

-5 

11 

-1 

2 

9 

-39 

-13 

11 

-6 

52 

17 

-31 

-2 

3 

-3 

4 

-12 

66 

521 

449 

-131 

14565 

-3 

4 

20 

-184 

10557 

2683 

-4 

5 

12 

199 

142 

1030 





1 



2 

-1 

1 

— 7 

17 

-2 

2 

9 

-47 

-16 

14 

-9 

51 

19 

-41 

-3 

3 

-3 

3 

-14 

85 

373 

533 

-81 

15367 

-4 

4 

17 

-229 

-8071 

2953 

-4 

4 

15 

233 

171 

1086 



Uz 







1 



2 

-3 

21 

-1 

3 

7 

-65 

-10 

17 

-5 

54 

13 

-47 

-2 

4 

-2 

4 

-9 

93 

443 

627 

-99 

19277 

-3 

5 

18 

-243 

-8883 

3509 

-3 

5 

10 

258 

135 

1299 



Trace of units, 
*(n), t(v~ l ) 



-1 
-2 

1 

-4 

2 

-5 

15 

-49 

4 

-7 

-25 

151 

39 

-42 

20 

-157 

-49 

119 

7 

-10 

8 

-11 

35 

-244 

-1337 

-1609 

311 

-49209 

10 

-13 

-55 

656 

27511 

-9145 

11 

-14 

-37 

-690 

-448 

-3415 



CO 
^4 



199 


66 


-4 


5 


-59 


21 


25 


20 


-46 


-41 


-46 


211 


70 


4 


5 


125 


25 


25 


20 


-46 


-45 


-50 


223 


74 


9 


2 


256 


28 


24 


22 


-47 


-50 


-52 


229 


76 


7 


4 


212 


28 


26 


22 


-49 


-50 


-54 


241 


80 


-6 


5 


-125 


25 


30 


25 


-56 


-50 


-55 


271 


90 


-10 


3 


-261 


27 


33 


30 


-64 


-57 


-60 


277 


92 


-9 


4 


-236 


28 


34 


30 


-65 


-58 


-62 


283 


94 


-11 


2 


-304 


28 


34 


32 


-67 


-60 


-62 


307 


102 


5 


6 


216 


36 


36 


30 


-67 


-66 


-72 


313 


104 


-12 


1 


-371 


31 


37 


36 


-74 


-67 


-68 


331 


110 





7 


49 


37 


40 


33 


-74 


-70 


-77 


337 


112 


-2 


7 


-25 


37 


41 


34 


-76 


-71 


-78 


349 


116 


12 


1 


517 


43 


37 


36 


-74 


-79 


-80 


367 


122 


-12 


3 


-435 


37 


44 


41 


-86 


-78 


-81 


373 


124 


4 


7 


221 


43 


44 


37 


-82 


-80 


-87 


379 


126 


-10 


5 


-365 


39 


46 


41 


-88 


-80 


-85 


397 


132 


11 


4 


544 


48 


44 


40 


-85 


-88 


-92 


409 


136 


10 


5 


515 


49 


46 


41 


-88 


-90 


-95 


421 


140 


6 


7 


343 


49 


49 


42 


-92 


-91 


-98 


433 


144 


-1 


8 


16 


48 


52 


44 


-97 


-92 


-100 


439 


146 


9 


6 


504 


52 


50 


44 


-95 


-96 


-102 


457 


152 


3 


8 


220 


52 


54 


46 


-101 


-98 


-106 


463 


154 


-8 


7 


343 


49 


56 


49 


-106 


-98 


-105 


487 


162 


8 


7 


505 


57 


56 


49 


-106 


-106 


-113 


499 


166 


-11 


6 


-536 


52 


60 


54 


-115 


-106 


-112 



33 over 
63 over 
35 over 
23 over 
27 over 
30 over 

34 over 
33 over 
30 disc 
39 over 

30 disc 
39 over 
57 over 

43 over 

35 over 
38 over 
60 disc 

31 disc 
31 disc 
50 disc 
49 over 
47 over 

44 over 
33 over 
56 over 



2f 

If 

If* 

If* 

If 

If* 

If 

If* 

If* 

4ff 

If 

If* 

7ff 

If* 

Is* 

Is* 

If 

If* 

If* 

If* 

If* 

If* 

Is* 

If* 
If* 



-45 

95 

17 

233 

-117735 

61785 

5673 

1 62609 

-107 

87 

-37 

6029 

-27 

27 

14945 

-38735 

-398895 

1233880499 

-5 

6 

-31 

29 

127 

5189 

-6 

7 

10967 

654389 

-863424 

-77129 

-269376 

-258046 

37 

-627 

5133 

495 

122 

1167 

-111055 

91329 

958722876509 

-2609272963 

-4710533 

2763599779 

-57726 

-49086 

10649 

2607 

4533791744955 

-3447237 



-35 

105 

12 

273 

-120801 

49591 

6737 

129785 

-86 

94 

-19 

6254 

-21 

291 

17561 

-39631 

-4265923 

1405925741 



-34 

32 

144 

5685 

-6 

6 

12600 

671994 

-924319 

-862345 

-304892 

-269967 

33 

-719 

5330 

418 

130 

1304 

-120569 

99383 

814793743515 

-2939073045 

-5036695 

3040707889 

-49333 

-51886 

11170 

2904 

5045136467209 

-3665999 



-39 

121 

16 

296 

-96959 

60217 

7111 

154103 

-93 

108 

-28 

7316 

-23 

337 

14607 

-46569 

-3500815 

1503532433 

-5 

7 

-28 

35 

116 

6266 

-5 

7 

10682 

771641 

-772262 

-923164 

-257481 

-305561 

35 

-741 

4511 

476 

109 

1381 

- 102055 

107897 

917779916345 

-3070187671 

-4281243 

3251216741 

-52147 

-57437 

9561 

3046 

4350456032733 

-4304363 



119 

-321 

-45 

-802 

335495 

-171593 

-19521 

-446497 

286 

-289 

84 

-19599 

71 

-905 

-47113 

124935 

11755695 

-4143338673 

16 

-19 

93 

-96 

-387 

-17140 

17 

-20 

-34249 

-2098024 

2560005 

2556806 

831749 

833574 

-105 

2087 

- 14974 

-1389 

-361 

-3852 

333679 

-298609 

-2691296536369 

8618533679 

14028471 

-9055524409 

159206 

158409 

-31380 

-8557 

13929384244897 

11217599 



11. Conclusion 

The EDVAC produced 176 pages of tabulations, too long to be reproduced here, including 
some 3,500 algebraic numbers with small norms, in less than 6 hours. With the precision 
provided for in the coding, only four cases ran to completion, but 25 of the 43 cases run produced 
at least one unit in the process (for Z=97, 13 units were produced). The remaining 18 cases 
provided enough material to permit the hand computation of at least one unit. (For /=499, the 
unit required 43 decimal places to check by hand). 

In every case where g=l the explicitly known fundamental unit (—(k—g-\-l)/2, 
— •(&—#+ 1)/2, —(k—g—l)/2 was among those produced by machine, as well as the (non- 
fundamental) unit (—1,-2,0), except when 1=37 and 313. 

In addition to units, mentioned earlier, the number of matrices produced after </> is shown 
in table 2, with the following symbols denoting the outcome of the computation (see section 7): 
"comp" meaning completed, "disc" meaning incomplete (by failure of sign discrimination), 
"over" meaning incomplete (by overflow in norm formation). 

It is clear, however, that to obtain the fullest benefit of this algorithm one would require 
considerably greater storage, both to perform the unmodified algorithm (guaranteeing a 
fundamental unit) and to permit sufficient decimal (or preferably integral) precision. 
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